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Abstract 

We consider splitting and stabilization of second-order solitons (2-soliton breathers) in a model 
based on the nonlinear Schrodinger equation (NLSE), which includes a small quintic term, and weak 
resonant nonlinearity management (NLM), i.e., time-periodic modulation of the cubic coefficient, 
at the frequency close to that of shape oscillations of the 2-soliton. The model applies to the light 
propagation in media with cubic-quintic optical nonlinearities and periodic alternation of linear 
loss and gain, and to BEC, with the self-focusing quintic term accounting for the weak deviation of 
the dynamics from one-dimensionality, while the NLM can be induced by means of the Feshbach 
resonance. We propose an explanation to the effect of the resonant splitting of the 2-soliton under 
the action of the NLM. Then, using systematic simulations and an analytical approach, we conclude 
that the weak quintic nonlinearity with the self-focusing sign stabilizes the 2-soliton, while the self- 
defocusing quintic nonlinearity accelerates its splitting. It is also shown that the quintic term 
with the self-defocusing/focusing sign makes the resonant response of the 2-soliton to the NLM 
essentially broader, in terms of the frequency. 

PACS numbers: 03.75.Lm; 42.65.Tg; 42.81.Dp; 05.45.Yv 
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A well-known consequence of the integrability of the one-dimensional nonlin- 
ear Schrodinger equation with the self-focusing cubic nonlinear term is that this 
equation gives rise to an infinite family of iV-th order soliton solutions, in the 
form of periodically oscillating breathers (the oscillation period is the same for 
all N > 2). These solutions may be interpreted as bound states of iV overlap- 
ping fundamental solitons with different amplitudes. Their binding energy being 
exactly zero, all iV-soliton complexes are unstable against perturbations of the 
initial conditions, which can split them into slowly separating constituent fun- 
damental pulses. In this work, we aim to study the stabilization/destabilization 
of the 2- (second-order) solitons by a weak additional quintic self-focusing/self- 
defocusing term. Various origins of such additional nonlinearity, of either sign, 
are known in Bose-Einstein condensates (BECs) and nonlinear optics, being 
within the reach of available experimental techniques. We argue that our pre- 
dictions of the change in the stability of the higher-order solitons, and of a 
possibility to implement the precise control of the soliton dynamics by means 
of the quintic nonlinear term, may be realized in experiments with optical and 
matter-wave solitons. The analysis is performed using a combination of direct 
simulations and analytical approximations. 



I. INTRODUCTION 



The control of soliton dynamics has been drawing a great deal of interest both as a 
fundamental problem and a topic with a vast spectrum of applications - in particular, to 
optics and matter waves It is well known that the nonlinear Schrodinger equation 
(NLSE), being an integrable one, supports an infinite sequence of exact higher-order soliton 
solutions, which may be understood as bound states of strongly overlapping fundamental 
solitons [2j . However, within the framework of the integrable equation, the binding energy 
of the exact multi-soliton complex is always equal to zero 2] , hence the higher-order solitons 
are unprotected (unstable) against perturbations of initial conditions, which can induce 
splitting into fundamental constituents. For example, the second- and third-order solitons 
(which are often briefly called 2-soliton and 3-soliton, respectively) readily split into sets of 
two and three fundamental solitons, with amplitude ratios 1 : 3 and 1:3:5. 
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Different approaches were proposed to stimulate the splitting and make it a real physical 
effect. One possibility is to introduce a specific nonlinear dissipation into the model by 
adding nonconservative terms to the NLSE, such as the one accounting for the intrapulse 
Raman scattering in optical fibers [3j. Other physically relevant settings are those with a 
periodic modulation of either the group- velocity-dispersion (GVD) or nonlinearity coefficient 
in the NLSE, which is known as the dispersion management (DM) and nonlinearity manage- 
ment (NLM), respectively jf|. The DM in optical fibers can be of both sign-alternating and 
sign-preserving types (in the former case, the GVD coefficient periodically changes between 
positive and negative, alias normal-GVD and anomalous-GVD, values). The format of the 
DM can be piecewise-constant, built as an alternation of fiber segments with different values 
of the GVD coefficient, or sinusoidal. It was predicted that the latter format, with a mild 
modulation amplitude that does not imply the change of the sign of the GVD, can induce 
the splitting of both fundamental [3] and higher-order solitons 0, 0]. Both effects have 
been experimentally demonstrated in Ref . (| , which made use of a specially fabricated fiber 
with the diameter subjected to the sinusoidal modulation along the fiber, thus inducing the 
modulation of the local GVD coefficient. 

Various effects of the NLM in the context of fiber-optic telecommunications were the- 
oretically studied in Refs. 7, and the integration of the NLM with the DM was 
considered in Ref. [h]]. In terms of Bose-Einstein condensates (BECs) in dilute atomic 
gases, which, in the mean-field approximation, is also described by the NLSE (called the 
Gross-Pitaevskii equation, in that context [111]), the NLM represents the application of the 
Feshbach-resonance technique to the BEC in the case when the resonance is induced by 
a variable (ac) magnetic field 12j. In the latter case, Ref. 14| has demonstrated that a 



small-amplitude variable part of the nonlinearity coefficient gives rise to resonant splitting 
of higher-order solitons into fundamental ones, provided that the NLM frequency is close to 
the frequency of free shape oscillations of the higher-order bound state. 

In this work we demonstrate that the addition of a very weak quintic nonlinearity dra- 
matically changes the stability of 2-solitons under the action of the NLM. Namely, the 
2-soliton's splitting time becomes significantly larger (smaller) in the presence of weak addi- 
tional attraction (repulsion), represented by a self-focusing (defocusing) quintic term. These 
predictions may be tested experimentally in nonlinear optics and BEC. A challenging possi- 
bility is to create effectively stable 2-solitons using a weak self-focusing (attractive) quintic 



3 



nonlinearity, and control their behavior by means of the weak NLM. While both the small 
quintic term and weak NLM represent perturbations that break the integrability of the un- 
derlying NLSE, the stabilization of the 2-soliton states under the combined action of both 
perturbations implies an intriguing possibility of their mutual compensation, leading to an 
effective extension of the quasi-integrable behavior of the solitons. 

In optics, the cubic-quintic (CQ) nonlinearity with different sign combinations of the two 
;erms were predicted 15] and observed [3] in aqueous colloids, as well as in dye solutions 
Q and, very recently, in thin ferroelectric ffl ms Q. The sanre nonhnearity was also 
predicted, via the cascading mechanism, in two-level media [19[. The self-defocusing quintic 
nonlinearity, accounted for by a proximity to the resonant two-photon absorption, was also 



observed in other optical materials 



20|. 



In the description of effectively one-dimensional BEC settings, a universal self- attractive 
quintic term in the respective Gross-Pit aevskii equation (as said above, the attractive quintic 
nonlinearity should facilitate the creation of stabilized 2-solitons) accounts for the deviation 
from the exact one-dimensionality, i.e., a finite transverse size of the corresponding trap for 



the atomic condensate 
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23(| . Besides that, a self-defocusing quintic term may take 



into regard three-body collisions in the BEC, provided that the related losses are negligible 



24]. 



As concerns the time-periodic modulation of the cubic coefficient, dealt with in this work, 
in optics it may naturally arise as a result of the periodic alternation of the linear loss and 
compensating gain (a well-known transformation removes the respective linear terms in the 
NLSE, mapping them into the effective NLM) p|. As mentioned above, the same modulation 
in BEC represents the action of the Feshbach resonance controlled by the ac magnetic field. 
Thus, both the CQ nonlinearity and NLM are generic features of numerous physical settings. 

The paper is organized as follows. Results obtained by means of systematic simulations, 
that demonstrate the stabilization/destabilization of the 2-soliton, subject to the action of 
the resonant or near-resonant NLM, under the action of the weak quintic term, that corre- 
sponds, respectively, to the self-attraction/repulsion, are reported in Section II. Analytical 
approximations, which make it possible to explain the underlying effect of the resonant 
splitting of higher-order solitons under the action of the weak NLM, and also the stabiliza- 
tion/destabilization of the 2-soliton, induced by the quintic term, are presented in Section 
III. These approximations are based on analysis of the system's energy in the presence of 
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the NLM and quintic term. Finally, Section IV concludes the paper. 



II. SPLITTING OF THE SECOND-ORDER SOLITONS 

A. Second-order soliton in nonlinear Schrodinger equation with cubic-quintic non- 
linearity 

In this work, we take the one-dimensional NLSE for wave function </> (x, t) in the usual 
scaled form, 

where g > is the coefficient accounting for the cubic self-attraction, which is set to be g = 1 
in the absence of the NLM. In physically relevant realizations of Eq. ([I]), the dimensionless 
quintic-interaction constant, which accounts for the higher-order self- attraction/repulsion 
in case of e < 0/e > 0, is a small parameter, |e| <C 1. Energy E (the Hamiltonian) , from 
which Eq. ([1]) can be derived as id<p/dt = 5E/5<f)*, where 5 /5(f)* stands for the variational 
derivative with respect to the complex-conjugate field |2|, is 



\-oo 

E 

oo 



dx 



2 



dx. (2) 



Exact fundamental-soliton solutions to Eq. ([I]) are known for either sign of e [23, [26|] . In 
the case of the integrable cubic NLSE, with e = 0, exact n-soliton solutions are generated 
by initial conditions 

(j)(x, t = 0)=nA sech(Ar), (3) 

with integer n > 1 j^ . For all n > 2, they are breathers whose shape oscillates with the 
frequency independent of n, 

uJ sh = AA 2 (4) 

(27r/o; s h is usually called the soliton period). The explicit solution for the 2-soliton is rela- 
tively simple: 

Six t) = \Ac iAHl2 C ° Sh ^ AX) + 3e4 ^ 2f C ° Sh ^ (5) 
VK ' ' cosh (4Ax) + 4 cosh (2 Ax) + 3 cos (4AH) ' 1 ; 

The energy of the 2-solution, taken as per Eq. ([2]) with g = 1 and e = 0, is 

(£o) 2 _ sol = -(28/3)A 3 (6) 




FIG. 1: Oscillations of the peak density of the second-order soliton in the absence of the quintic 
nonlinearity, e = 0, (a), and in the presence of a weak attractive quintic nonlinearity, with e = 
-1.6 x 1(T 3 , (b). 

(as said above, it is identical to the sum of the energies of two far separated fundamental 
solitons with amplitudes 3A and A into which the 2-soliton complex may split). Figure Q^a) 
displays oscillations of the peak density \<p{x = 0,t)| 2 of the 2-soliton, generated by initial 
condition (J3]) with n = 2. 

Numerical simulations of Eq. ([I]) are performed using the split-step Fourier spectral 
method with 4096 Fourier modes and absorptive boundary conditions. Figure [T](b) displays 
504 oscillations of the peak density of the numerical solution with a weak attractive quintic 
nonlinearity, corresponding to e = —1.6 x 10~ 3 . Despite the absence of exact solutions for 
higher-order solitons at e 7^ 0, the oscillations persist indefinitely long without any tangible 
decay or distortion, the only difference from the exact solution presented in panel (a) for 
e = being a small increase in the oscillation frequency, as a result of the weak additional 
self- attraction. In view of the apparent robustness of this solution at small values of |e|, we 
will continue to refer to it as the second-order soliton. 
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Position 



FIG. 2: A typical example of the splitting of the second-order soliton under the action of a weak 
resonant modulation of the coefficient in front of the cubic nonlinearity, in the absence of the 
quintic term. 

B. Resonant splitting of the second-order soliton 

We introduce the NLM, in the form of the time-periodic (ac) modulation of the nonlin- 
earity, by setting 

g(t) = l + bsm(ut), (7) 

in Eq. Qj, where amplitude b of the perturbation is small, |6| <C 1, and the modulation 
frequency is kept in resonance with the shape oscillations, to = u^. Figure [2] shows a 
typical example of the evolution of the wave function subject to the action of the resonant 
perturbation for e = (with the cubic nonlinearity only) and b = 5 x 10~ 3 . The 2-soliton 
splits into fundamental solitons with amplitudes related expected from the exact 

solution available at b = 0. The respective velocity ratio of the splinters is 3 : 1, in agreement 
with the prediction based on the momentum conservation (it follows from the fact that the 
effective masses of the two splinters are in the same ratio as their amplitudes [2J, i.e., 1 : 3, 



and the total momentum must remain equal to zero [l4|). 

The simulations demonstrate that even a very weak resonant perturbation splits the 
second-order soliton, but the onset of the splitting requires time which increases with the 
decrease of the strength of the ac modulation, b. The smallest perturbation amplitude that 
we applied in the simulations was b = 5 x 10~ 5 . This limitation was determined by the 
accumulation of numerical errors and available computational time; however, the results 
strongly suggest that there is no cutoff for the resonant splitting, as concerns the smallness 
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FIG. 3: A detailed picture of the decay of density oscillations in the course of the splitting of a 
second-order soliton into its two fundamental constituents. Oscillations of the highest peak in the 
two-soliton waveform are shown during and after the splitting. 

of the perturbation amplitude. 

For the quantitative analysis of the splitting, it is necessary to exactly define the splitting 
time, T. Figure [3] shows a generic example of the gradual decay of density oscillations of 
a second-order soliton in the course of its splitting. Following this picture, we define T as 
the time elapsed from the start of the action of the NLM until the amplitude of the density 
oscillations drops to 10% of its initial value in the 2-soliton. By that time, the second-order 
soliton turns into well-separated constituents. In fact, adopting another technical definition 
of the splitting time produces a little effect on the final results. 

The chain of open squares in Fig. H] shows the splitting time as a function of the pertur- 
bation strength, b (in the case of e = 0), revealing the divergence in the limit of b — > 0. We 
fit this dependence to a simple power-law expression, 

T = a ■ V + c (8) 

where a, c and p < are constant parameters and b is taken in percents. As can be seen in 
the figure, the fitting formula captures the behavior of the numerical data well. For e = 0, the 
fitting parameters are p = -0.462±0.022 and a = (5.24±0.59) x 10 3 , c = (1.44±0.78) x 10 3 . 

Further, black triangles and circles in Fig. H]show the splitting time as a function of the 
NLM strength in the presence of the weak quintic nonlinearity of either sign, viz., for 

e = -9.1 x 1(T 4 and e = 3.03 x 10~ 4 , (9) 
8 
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FIG. 4: The splitting time of the 2-soliton, as produced by the simulations, versus the ac-drive's 
amplitude, b. Open squares pertain to the case of the cubic nonlinearity only (e = 0), black 
triangles are for e = —9.1 x 10~ 4 , and black circles are for e = 3.03 x 10~ 4 . Solid lines show 
the fit provided by Eq. ([8]). For positive e, the fitting formula (shown by the dashed line in this 
region) does not describe the numerical data at smallest values of b, which indicates saturation for 
extremely weak perturbation strengths. 

respectively. Even this very weak extra nonlinearity affects the splitting time dramatically: 
the additional quintic attraction (repulsion) delays (accelerates) the splitting, thus effectively 
stabilizing (destabilizing) the 2-soliton. In particular, the quintic self-focusing term with the 
amplitude as small as e = —9.1 x 10~ 4 (black triangles in Fig. H]) is enough to increase the 
spitting time by a factor of ~ 2.5 when the ac-drive's strength is b = 1 x 10~ 3 , as compared 
to the case of e = 0. 

The fit parameters for the data pertaining to the extra quintic attraction or repulsion, 
with the values of e as in Eq. flSj), are, respectively, 

p= -1.93 ±0.16, a = (4.2 ±1.5) x 10 2 , c = (7.25 ± 0.53) x 10 3 , (10) 

p = -0.372 ± 0.007, a = (5.84 ± 0.09) x 10 3 , c = 0. (11) 

It is worth noting the differences in parameter c for both cases. For e < 0, positive c in the 
fitting set ([TO]) means that, even for a strong perturbation, a final waiting time is required 
to observe the splitting of the 2-soliton, which is another manifestation of its stabilization 
by the quintic self-focusing term. On the contrary, for e > 0, the best fit actually required to 



choose c < 0-in the parameter region where formula (jHJ) with c < produces T > 0. Setting 
c = in the fitting set (JTTJ) means that the simulations demonstrate that the 2-soliton starts 
splitting instantaneously under the action of the strong resonant perturbation. 

Actually, for large perturbation amplitudes (b > 0.02), the splitting produces fundamental 
solitons with the amplitude ratio different from 1 : 3, which may be explained by effects 
induced by the relatively strong perturbation on the constituents of the 2-soliton in the 
course of the splitting. In the case of the self-defocusing quintic nonlinearity, e > 0, the 
splitting time shows saturation for extremely weak perturbations (see black circles in Fig. 
HJ), i.e., the splitting time ceases to grow with further weakening of the perturbation. A 
plausible explanation to this feature, which demonstrates the fragility of the second-order 
soliton in this situation, is the fact that splitting is spontaneously initiated by the numerical 
noise. We also note that the quintic nonlinearity slightly changes the frequency of the shape 
oscillations of the 2-soliton, see Fig. [IJ The modulation frequency was modified, accordingly, 
in the simulations, to maintain the resonance condition for all cases included in Fig. HI 

Figure [5] summarizes absolute values of exponent \p\, i.e., the best-fit parameter in Eq. 
(jHJ), as a function of strength e of the quintic nonlinearity. The error bars in the figure 
represent errors of the fit to the power-law approximation defined as per Eq. (jSj). For 
e < 0, \p\ rapidly increases with — e, indicating that the second-order soliton develops strong 
resistance against the resonant splitting, with the enhancement of the quintic self-focusing. 
On the other hand, \p\ drops for e > 0, showing the destabilization of the second-order 
soliton under the action of the quintic self-defocusing. 



C. The near-resonance response 



In Ref. 14j it was shown that the splitting of the 2-soliton (in the case of e = 0) could also 
be caused by the temporal modulation of the coefficient in front of the cubic term with the 
frequency slightly different from resonant value (j3J). Here, we aim to confirm this behavior 
for case of the cubic NLSE and extend the analysis to the CQ model, with e ^ (in the 
latter case, the resonant frequency should be first slightly adjusted, as mentioned above). 
For this purpose, we have performed the simulations by varying modulation strength b at 
fixed values of the driving frequency and fitting the so observed splitting time (T) to the 
power-law function, defined as per Eq. (jHJ). 
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FIG. 5: The power-law fit parameter, \p\, from Eq. (jHJ), as a function of the strength of the quintic 
nonlinearity, e. The error bars represent the error of the fit to formula (jHJ). 

Figure M displays the absolute value of the fitting power, \p\, versus the modulation 
frequency, for four different values of the quint ic-nonlinearity coefficient, e = —1.6 x 10 -3 
(open rhombuses), e = —7.6 x 10~ 4 (black triangles), e = (open squares), and e = 7.6 x 10~ 4 
(black circles). As before, the error bars in the figure represent errors of the fit. Two 
noteworthy features are clearly seen in the figure. First, small and large values of \p\ are 
obtained for the repulsive and attractive quintic nonlinearity, respectively, in accordance 
with what was reported above for the case of the exact resonance. Second, in most cases the 
presence of the quintic term of either sign gives rise to broadening of the resonant response, 
in comparison with the case of the pure cubic nonlinearity (represented by open squares 
in Fig. [6) the broadening is weak in the case of the weakest quintic self-focusing, which 
corresponds to e = —7.6 x 10 -4 , which is represented by the chain of black triangles). In the 
case of the repulsive quintic nonlinearity, e > 0, the dependence shown by the circles in Fig. 
[6]is nearly flat, i.e., the second-order soliton readily splits even at a relatively large detuning 
from the resonance. On the contrary, the stabilization of the 2-soliton by the self-focusing 
quintic nonlinearity is seen to be robust also under the off-resonance conditions. 
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FIG. 6: The power-law fit parameter, \p\, as a function of the modulation frequency. Open squares, 
black triangles, open rhombes, and black circles are, respectively, for the pure cubic nonlinearity 
(e = 0), e = —7.6 x 10~ 4 , e = —1.6 x 10~ 3 (the self-focusing quintic term), and e = 7.6 x 10~ 4 (the 
self-defocusing one). The error bars represent errors of the fit to the power-law expression defined 
as per Eq. (jSj). 

III. ANALYTICAL ESTIMATES 

A qualitative explanation to some numerical findings reported above can be provided by 
an analytical consideration of the model based on Eq. ([1]). First, it is possible to explain the 
underlying effect of the resonant splitting of the 2-soliton. Indeed, as mentioned above, the 
binding energy of the higher-order solitons is exactly zero in the integrable NLSE, therefore 
the splitting may be explained by the fact that the resonant temporal modulation pumps 
energy into the bound 2-soliton state, causing its splitting into constituents, which carry 
away the excess energy, in the form of their kinetic energies. In the presence of the small 
modulation term in the cubic coefficient given by Eq. ([7j), whose frequency is set to coincide 
with the resonant value (j3J, the exact evolution equation for energy Eq of the unperturbed 
NLSE, i.e., Eq. (j2J) with g = 1 and e = 0, can be derived in the following form: 




In the lowest approximation, one can substitute the unperturbed 2-soliton solution, as given 
by Eq. (jSJ), into the right-hand side of Eq. (fl2"|) . Under the condition that the shape 
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oscillations of the 2-soliton are synchronized with the temporal modulation in Eq. (J7J) 
(analysis of of numerical results confirms this conjecture, in the case of the slowly developing 
splitting), the time averaging of Eq. (TT21 yields an effective energy-pump rate, 

I dE \ 128 



, , , -C b A % « 23.2 A%, 
\ at I n 

where the constant is given by the following integral expression, 



(13) 



a 



2- 



dr sin r 



<i£ Im 



cosh (3£) + 3e lT cosh £ 
cosh (4f ) + 4 cosh (2f ) + 3 cos r 



9 cosh (3£ ) + 3e lT cosh £ 



n 2 



0.57. 



(14) 



<9£ cosh (4f ) + 4 cosh (2f ) + 3 cos r 

This explanation of the gradual onset of the splitting of the 2-soliton through the pumping 
of the energy into it by the resonant NLM was not considered in Ref. [14], which was dealing 
with the (near-)resonant splitting in the cubic NLSE (with e = 0). 

The stabilization/destabilization of the 2-soliton by the self-focusing/defocusing quintic 
nonlinearity may also be explained by means of the consideration of the energy. In expression 
([2]), the term corresponding to the quintic nonlinearity in Eq. (CQ) yields the following 
expression for the additional energy, obtained by the substitution of the unperturbed 2- 
soliton solution ([5]) and averaging over the period of its shape oscillations: 



(A3 



2— soliton 



) = (2 n a/3vr) A 5 



39.5 A 5 e, 



(15) 



with constant 

C e = dr 
Jo 



{[cosh (3£) + 3 (cos r) cosh £] 2 + (sin 2 rj cosh 2 £ j 
d£, ; : — ; — : : — ; — — ~ 0.18. 



[cosh (4f) + 4 cosh (2f) + 3 cos rf 
Finally, the binding energy of the 2-soliton induced by the quintic term can be defined as 
the difference between the respective extra energies of the pair of free fundamental solitons, 
with amplitudes 3A and A, into which the 2-soliton splits, AE 3 a+a — (3 5 + 1) (16/45)eA 5 rs 
86.8 eA 5 , and energy ( fl~5l) of the unsplit bound state: 



AE = AE 3A+A - AE 2 _ sohton ~ 47.3 eA 5 . 



(16) 

The negativeness and positiveness of expression (f!6l) for e < and e > 0, respectively, 
explains the stabilization/destabilization of the 2-soliton by the self-focusing/defocusing 
quintic nonlinearity. 
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According to Eq. ([13]) . under the action of the resonant NLM the energy of the resonantly 
driven 2-soliton grows, on the average, linearly in time: 

E (t) « E (0) + 23.2 A 5 bt, (17) 

where the initial value, Eq(0), is energy ([6]) of the unperturbed 2-soliton. If all the energy 
pumped by the NLM drive would be accumulated in the form of the "splitting potential" , 
then, in the case of e < 0, the onset of the splitting might be expected at the moment of 
time when this stored "potential" becomes equal to the binding energy, which would give 
the splitting time as T = AE/ (dE /dt). In its literal form, this result implies p = — 1 in 
Eq. (jSj). The deviation of actual values of \p\ from 1, as seen in Fig. [5j is explained by the 
fact that, within the framework of the present analysis, it is not known how the pumped 
energy is divided between the actual accumulation of the splitting potential and absorption 
into a change of the unperturbed 2-soliton's energy - see Eq. (jHJ) - due to a possible small 
variation of A. While an exact prediction of \p\ as a function of e seems to be too difficult 
for the fully analytical treatment, the fact that, according to Fig. [5l \p\ is essentially larger 
than 1 at — e > 10~ 3 suggests that, in this range of values of e, nearly all the pumped energy 
is absorbed into the increase of the 2-soliton's energy. Combining Eqs. (jUJ) and (fTTj) . one 
can conclude that, as long as the resulting deviation of amplitude A from its initial value, 
Aq, remains small, the amplitude varies in time as A(t) ~ Aq — 0.8Aq6 • t. This variation 
leads to a detuning of the resonant driving, according to Eq. (TjJ, but, on the other hand, 
Fig. [6] demonstrates that the detuning does not produce an essential effect for — e > 10~ 3 . 



IV. CONCLUSIONS 



In this work we have considered the influence of the weak quintic nonlinearity on the 
stability and splitting of second-order solitons (alias 2-solitons) in the NLSE-based model, 
under the action of the weak resonant NLM (nonlinearity management), i.e., periodic time 
modulation of the coefficient in front of the cubic term, with the frequency equal or close 
to the frequency of the free shape oscillations of the 2-soliton. The model applies to the 
propagation of light in CQ optical media, taking into regard the periodic action of the linear 
loss and compensating gain. The same model finds a natural application to BEC, where 
the self-focusing quintic term accounts for the effect of the residual three-dimensionality 
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in the effective one-dimensional approximation, while the NLM may be induced by the 
Feshbach resonance controlled by an ac magnetic field. By means of direct simulations and an 
approximate analytical considerations, we have demonstrated that the additional weak self- 
focusing quintic nonlinearity stabilizes the 2-soliton, while the self-defocusing nonlinearity of 
the same type makes it fragile and accelerates its splitting. We have confirmed the resonant 
character of the splitting of the 2-soliton under the action of the NLM, and proposed an 
explanation to this effect, based on the consideration of the rate at which the energy is 
pumped into the bound state by the resonantly tuned ac drive. We have also studied 
the resonant NLM-induced splitting of the 2-soliton in the presence of the weak quintic 
nonlinearity. Depending on its sign, the self-defocusing/focusing higher-order nonlinearity 
gives rise to conspicuous broadening/sharpening of this resonant response. The results of the 
numerical and analytical considerations reported in this paper for 2-solitons can be readily 
extended to higher-order solitons. 
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